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Abstract. In this paper, we prove the global existence and singularity formation for a 
wave system from modelling nematic liquid crystals in one space dimension. In our model, 
although the viscous damping term is included, the solution with smooth initial data still 
has gradient blowup in general, even when the initial energy is arbitrarily small. 
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CN) 1. Introduction 

The mean orientation of the long molecules in a nematic liquid crystal is described by a 
director field of unit vectors, n G S 2 , and the propagation of the orientation waves in the 
director field could be modelled by below Euler-Largrangian equations derived from the least 
^ action principle [121 Q] > 

1.1 n u + fxn t + y ^ = An, n • n = 1, 

i ^ i on 

where the well-known Oseen-Franck potential energy density W is given by 

(1.2) W (n, Vn) = |a(V • n) 2 + \f3 (n • V x n) 2 + ^|nx(Vx n)| 2 . 

^ The positive constants a, /3, and 7 are elastic constants of the liquid crystal, corresponding 

to splay, twist, and bend, respectively. The viscous coefficient fj, is a non-negative constant 
and the Lagrangian multiplier A is determined by the constraint n • n = 1. 

There are many studies on the constrained elliptic system of equations for n derived 
through variational principles from the Oseen-Franck potential (1.2), and on the parabolic 
flow associated with it, see [21 O [[5] . 

The global weak existence and singularity formation for the Cauchy problem of the extreme 
^ case of (1.1) with \i = in one space dimension (1-d) has also been extensively studied 



O [H [16, [T71 [181 (HI SJ- This hyperbolic system describes the model with viscous effects 
neglected, for which an example with smooth initial data and singularity formation (gradient 
blowup) in finite time has been provided in [8 J . The lack of regularity makes us only be able 
to consider the existence of weak solutions instead of classical solutions. The existing global 
existence results in [31 [181 E3 S] proved by the method of energy-dependent coordinates will 
be introduced in Section [2l 



However, the well-posedness and regularity of the solution for the complete system (1.1) 
with fj, > are still wide open. 
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In this paper, we consider the 1-d system of (1.1) with a = (3 and fi > 0, where the 
derivation will be delayed to Section [2] 

d tt ni + fid t ni - d w {d t {nx)d x n 1 ) = (-\n t \ 2 + (2c 2 - 7)|n x .| 2 )m, 
dttn 2 + ^d t n 2 - d x (c 2 (m)d x n2) = (-\n t \ 2 + (2c 2 - a)\n x \ 2 )n 2 , 
dttn-z + fid t n 3 - d x {c 2 {ni)d x n 3 ) = ( — I n f j 2 + (2c 2 - a)\n x \ 2 )n 3 , 

c 2 (ni) ^ a + (7 — a)n\ 
and 

(1.4) n = n(t, x) = {n\(t,x),n 2 {t, x),n 3 (t,x)) and |n| = 1, where (t,i)el + xR. 

In this paper, we first provide an example with smooth (C 1 ) initial data and singularity 
formation in finite time. Then we prove the global weak existence of the Cauchy problem for 
dl.3|) with initial data 



(1.5) m(Q, x) = m £ H 1 , (nj)t(0, x) = nn £ L 2 , i = 1 ~ 3. 

Although when \x = the singularity formation and global existence of Cauchy problem 
for (1.3) has been systematically studied in [8] and [3] respectively, the extension of these 
results to the case with fi > in this paper is important and non-trival. First, the model 
with fi = only describes the extreme case in physics which is much less happening than the 
one described by the model with fi > 0. Secondly, the comparison between our results and 
the results when [i = indicates that the regularity of the solutions is basically not impacted 
by adding or removing the viscous term /Mt even when fi is large, which is to some extend 
unexpected. Furthermore, now most existing existence and regularity results for the models 
with fi = can be expected to extend to the models with \i > 0. 

We first consider the regularity of solutions of (1.3). A nature question is whether the 
appearance of the viscous term /xrit with /i > can prevent the singularity formation (gradient 
blowup), which is known existing when fi = 0, or not. However, the singularity formation 
result in this paper gives a negative answer to this question: For any [i > 0, the solution 
generally includes gradient blowup. In fact, for both the singularity formation example in [8] 
(// = 0) and the one in this paper (// > 0), the initial energy can be arbitrarily small. As 
a consequence of this surprising result, we need to consider the weak solution instead of the 
classical solution for system ( |1.3[ ). 

For simplicity, in order to get a gradient blowup example, we only consider the solution of 
(1.3) with structure n = (cos u(t, x), sinu(t, x), 0) (planar deformation). Clearly 713 is always 
in the smooth solution if it is vanishing initially. It is easy to get that the system (1.3) is 
equivalent to 

(1.6) utt + f-iut - c(u)(c(u)u x ) x = 0, 
when sinu 7^ 0, where we still use c to denote the wave speed and 

(1.7) c 2 (u) = 7 cos 2 ti + a sin 2 n. 

It is easy to see that there exist positive constants Cl, Cjj and Cd, such that 

(1.8) C L <c(u)<Cu, \c'(u)\<C D . 



Theorem 1.1. We consider the Cauchy problem of (1.6) with initial data satisfying 

U(0,X) = Uq + E(j)(^) + E 2 T](£^X), 

(1.9) u t (0,x) = (-c{u(0,x)) +e) u x (0,x), 
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where uq is a constant satisfying 

c'(u ) > 0, 

and two functions 4>{a) and 77(a) are in C 1 (M) and satisfy conditions 

(1.10) 0(a), 77(a) = when a £(—1,1); <j)'(a) < when a E [0, 1), 

and 

8/iCf/ 



(1.11) 



- 0'(O) > 



c'(no) C L 

For any fj,> 0, we can choose e > sufficiently small, such that the C 1 solution u(t, x) forms 
singularity in finite time. 



Note, in the solutions we constructed, sinu 7^ before blowup. Hence, for system (1.3), 



for any /i > 0, we can find examples with C 1 smooth initial data and singularity formation 
in finite time. From Theorem |1.2[ we know that the gradient blowup we found is not a 
discontinuity (shock wave). 

Then we consider the global existence of weak solution for the initial value problem of 
(1.3)~(1.5). Here, the unit vector n(t,x) with (t,x) E M + x 1 is a weak solution to the 
Cauchy problem (1.3)~(1.5) if it satisfies that 

1. The equations (1.3) hold in distributional sense for test functions (j) £ C*(R X R). 



The initial data satisfy (1.5) in the pointwise sense, and their temporal derivatives 
hold in L^ oc for p E [1,2). 
2. m(t,x) is locally Holder continuous with exponent 1/2, for i = 1 ~ 3. The map t 1— > 
(rai, 77,2, 77-3)(t, •) is continuously differ entiable with values in L?, for all 1 < p < 2. 
And it is Lipschitz continuous under the I? norm, i.e. 

(1.12) \\ni(t, •) - rii(s, -)\\ L2 < L \t - s\, i = 1 ~ 3, 

for all t,s e M+. 
The main well-posedness results in this paper are 



Theorem 1.2. The initial value problem (1.3)~(1.5) has a global weak solution n(t,x) for 
all (t, x) E M + x R, where n(t, x) satisfies that the energy 



(1.13) 



£(t) 



defl 



\n t \ 2 + ^(ni)!^ 



is less than or equal to £(0) for any t > 0. 



Theorem 1.3. For the Cauchy problem (1.3)~(1.4), let a sequence of initial data satisfy 



i n io)x 



till 



n io)x\\ L 2 ~^ 0, 

no uniformly on compact sets, as k 
subsets of the (t,x) -plane with t > 0. 



and rig 



00. 



i2 ^0, 
Then n k 



= 1 ~ 3, 

n uniformly on bounded 



In order to prove these theorems, by the method of energy-dependent coordinates, first 
used in papers [3] and related Camassa-Holm equation, we dilate the singularity then find 
the energy dissipated solution. In fact, the energy is dissipated when a < 0. When \i = 0, 
the energy is conserved in some sense [3]. 



Similar existence and continuous dependence results apply to system (1.6), which is a 



special example of (1.3) 



The paper is organized as follows. In Section 2, we derive the system (1.3) and energy 
equation, then introduce the existing results for the extreme case /j = 0. In Section 3, we 
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prove the singularity formation result: Theorem In Section 4, we construct a set of 
semi-linear equations in the energy-dependent variables by studying the smooth solution of 
(1.3). Then, in Section 5, we prove the existence of the solutions for the semi-linear system. 
In Sections 6~8, we transform the solutions back to the original system and prove Theorems 
OMOl 

2. The derivation of system 
In this section, we derive ( fL3l ) from ([IT]). In one space dimension, ( |1.1[ ) is 

(2.1) durii + fidtUi + d ni W(n, d x n) - d x [d dxni W(n, d x n)] = An;, for i = 1, 2, 3. 

Using |n| = 1, we get by multiplying rij to ( |2.1[ ) and summing up i from 1 to 3 

3 

(2.2) X = Y^[-\dtni\ 2 + md ni W(n, d x n) - md x [d dxrii W(n, d x ii)] }. 

i=l 



It is easy to calculate that (2.1) has the following energy equations: 

1 3 
(2.3) d t [-\n t \ 2 + W{n,d x n)] - d x [J2 dtmd dxni W (n, d x n)] = -^\n t \ 2 < 0. 



i=l 



In 1-d case, the Oseen-Franck potential energy density (1.2) is 

(2.4) W(n, d x n) = ^{d x n x f + ^ [{d x n 2 f + (d x n 3 ) 2 } + i( 7 - fJ)n 2 \d x n\ 2 , 
from which and (2.2), we infer 

(2.5) A = -|n t | 2 + [p + 2( 7 - ( 5)n?) \d x u\ 2 + (/3 - a)mSgm. 



Then by \2.l\, \2A\ and (|2.5|), we have 
(2.6) 



' d tt ni - dx [(%(ni)d x ni] = { — | n t | 2 + (2c| -7)|n :r | 2 + 2 (a - /3)(9 :c m) 2 }rai 

3 t in 2 + M n 2 - ^ [02(71-1)^77,2] = { — I n t 1 2 + (2c| - /3)|n x | 2 + (/3 - a)n 1 d xx n 1 }n 2 
> ^ttn 3 + /i3i«3 - ^[02(711)^713] = { — I n t 1 2 + (2c2 - /S))^! 2 + (/3 - a)nid xx ni}n 3 



with 



c 2 (ni) = a + (7 — a)n\ and c|(ni) = /3 + (7 — /3)n 2 . 



In particular, taking a = /3 in (2.6), we get (1.3) and 

(2.7) c 2 (ni) = c\ (ni) = c|(ni) = a + (7 - a)n 2 , 

and it is easy to check that there are positive numbers Cl < Cjj such that 

(2.8) < Cl < c(ni) < Cy < 00; |c'(ni)| < Cat, for any |m| < 1. 



The energy equation (2.3) is 

(2.9) \d t [|n 4 | 2 + c 2 (ni)| n;c | 2 ] - d x [c 2 (ni)n t • nj = -/x|n t | 2 < 0, 



where the energy density W in (2.4) is 

(2.10) W(n,d x n) = ^c 2 (n 1 )\d x n\ 2 . 

Finally, we introduce the existing existence results for weak solutions of 1-d systems of 



(1.1) with /x = 0, proved by the method of energy-dependent coordinates. This method was 
first applied to equation (1.6) with fi = for an existence proof by [3j. When n is arbitrary 
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in S 2 and fi = 0, under the a priori assumption that n(x,t) is uniformly away from (1,0,0), 



the existences of 1-d solutions for (1.3) and (2.6) with (3 < a have been provided by |18| and 



[19] , respectively. In a recent paper [3], the existence of 1-d solution for system (1.3) and 
(j, = has been established, without any a priori assumption. 



3. Singularity formation: Proof of Theorem 11.1 
We introduce Riemann variables R and S for equation ( |1.6[ ): 

R = f u t + c(u) u x , S = f ut - c(u) u x . 



(3.1) 

By (IX .6b , the C 1 solution satisfies 



(3.2) 
(3.3) 



Rt — cR x 
S t + cS x 



^(^-HVf(« + 5), 



where the left hand sides are two directional derivatives along characteristics. 

Before we give the proof, we first explain the basic idea. The singularity is essentially 
caused by the quadratic increase of the Riemann variables. In order to find it out, we need 
to get a decoupled Riccati type inequality on some characteristics. First, we show that the 
energy is in 0(e) for our example. Using this estimate, we can make sure that c!(u) > \d (uq) 
before t = a/e for some positive constant a, when e > is sufficiently small. This means that 
the coefficients of the quadratic terms in (3.2) (3.3) have fixed sign. Then the key step of this 



proof is to prove the domain: > R > —M^es, S > is invariant, when (t,x) is taken on 
some characteristic trapezoid. Collecting all these informations, we can find the decoupled 
Riccati type inequality then prove the blowup. 
We will prove the Theorem |1.1| in several steps. 

1. It is easy to get from the initial condition that 

,x 



(3.4) 

(3.5) 
(3.6) 



«a,(0,x) = 4>'(-) +e3 r]'(e3x), 

R(0,x) = eu x (0,x), 

5(0, x) = (-2c(w(0,s)) +e) u x {0,x). 



Choosing e < Cl and < £ < 1, we have, when x G [0, e a ) 
(3.7) u x (0,x)<0, R(0,x)<0, S(0,x)>0, 

because <j)'(a) and r/(a) are both negative when a G [0, 1). 



2. By (3.2) and (3.3), we have the energy equation 

(3.8) (R 2 + S 2 ) t + (c(5 2 - R 2 )) x = -n(S + R) 2 < 0. 



This equation agrees with the energy equations (2.3) and (2.9). We define the energy function 
E as 

/oo 1 Z" 00 

u 2 (x,t) + c 2 (u(x,t))u 2 x (x,t)dx = - / R 2 (t,x) + S 2 (t,x)dx. 
-oo * J~oo 
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Here we use E to denote energy for equation (1.6) instead of £ for (1.3) to avoid confusion 



Integrate ( |3.8| ) with respect to x, then we have 
E(t) < E(0) 

= - / R 2 (0,x) + S 2 (0,x)dx 

2 J-oo 

= \ r [{-2c(u{Q,x))+ef + e 2 ]{<!,' { X )+ sir,' (el x)fdx 

2 J-oo e 

< Me, 

for some positive constant M, where we use < e <C 1. 



(Xn.to) 




Figure 1. Characteristic triangle 



3. We consider any characteristic triangle in Figure [T] with the characteristic boundaries t± 
denoted by 

dt±(x) , 1 



d:X 



c(u) 



Integrating (3.8) on this characteristic triangle, and by the divergence theorem, we have 

rx 2 



(3.10) 



R 2 (t + (x),x)dx + S 2 (t^(x),x)dx < - 



.I'D 



%2 

[R 2 (0,x) + S 2 (0,x)] dx. 



4. Then we control the sign of c'(u) by choosing e small enough. Since equation ( 1.6 ) has finite 



propagation speed, clearly seen from (3.10), we get that u = uq in the region x < — e a — Cut 
and + C^. 



\u(t, X) — Uq\ 



u x (t, x) dx\ 



£ 3 +C V t 



-e Ti-Cijt 



\u x (t, x)\ dx 



(3.11) 



< 



< \\u x \\ L 2 \J2£-3 +2 Cut 



< 



C L 



2es +2C v te, 
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where we use the bound on E(t) in part 2. Hence, we can find a small positive number a 
independent of e then another small positive number eq, such that, 



(3.12) 
if 



c'(u(x,t)) > >0, 



< t < a/e, < e < e . 



We only consider the problem in this time interval, in which (3.12) is satisfied. 

5. We next prove an a priori estimate under the assumption that R < and S > in the 
characteristic triangle in Figure [lj with X2 — x\ < £~s. 
Under this assumption, we have 



Rt — cR x > 



S. 



Integrating it along the backward characteristic by (3.10), we have 

'to . rto 



R(t ,x ) 



> 



> 



> 



C D 
~±C L J 

c D 

Cp 
'4CP L 

-M 2 e\ 



S 2 (t,x-(t))dt 
X2 S 2 (t_(x),x) 

x C(t-(x),x) 







dx 



S(t,x-(t))dt + R(0,X!) 
X2 S(t_{x),x) 



c(t-(x), x) 



dx + R(Q,xi) 



E(0) 



2C L 



^/E^)^\x~^x~o~\ + R(0, xi) 



for some positive constant Mi. Choosing e < ( jj^cw ) 6 ' we nave 

^>0 



c '( u ) P 2 

; — r-ft 

4c(n) 2 



because R < under the a priori assumption we assumed. So 

- M ' 

4c(-u) " 



S t + > ^H^ 2 



-S. 



6. Then we use this a priori estimate to prove that R < and S > in the closed region 
bounded by the forward characteristic starting from the origin on (t, x)-plane, t = a/e, t = 
and the backward characteristic starting form the point (0, e _ s — s 2 ) on (t,x)-plane. Here, 

_ 2 r> _ 2 

we choose e 3 — e instead of e 3 to exclude the backward characteristic starting from the 

_2 

point (0, e 3). 

We prove it by contradiction. Assume that R = or S = at some point in the region we 
consider. By (3.7), R(0,x) < and 5(0, x ) > when x G [0, £~a], so we can find the lowest 
line t = t\ > such that R = or S = on some point of this line, while R < and S > 
when < t < tj.. So we can find a closed characteristic triangle with R = or S = on the 
upper vertex while < and S > elsewhere. Then we derive a contradiction by proving 
that R and S cannot be zero on the upper vertex. 

Since S > in the characteristic triangle, 

c '( n ) D 2 A 4 D 

Rf — cR r < - — ; — r-R — — XL. 

~ 4c(u) 2 



s 
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Then by standard ODE comparison theorem (see [IT] for reference), (1.8) and R < except 
the upper vertex, we have R < on the vertex. By S > and R < in the characteristic 
triangle, the a priori estimate in the previous part shows that 



S t + cS x >p^S*-»S. 
~ 4c(«) 2 



Then by ODE comparison theorem, (1.8) and S > except the upper vertex, we have S > 



on the vertex. Hence, we get a contradiction. 



7. Finally, we prove the singularity formation by considering the forward characteristic start- 
ing from the origin. We have already proved that along this characteristic, before t= — and 

- 2 9 



2C V 



(3.13) 



St + cS x > 4V - > C ^s" -^s> *M s 2 - Us. 

- 4c(«) 2 ~ 8C7{/ 2 ~ 16Q; 2 



When £ < C^, (3.6) and the initial conditions (1.10) and ( |1.11|) give that 
(3.14) 

which shows 



CW(0) > S(0,0) = -c(»(0,0))*'(0) > ^ 



^(0,0) > 25(0,0). 



By the ODE comparison theorem, on the characteristic, (3.13) gives 



c'{u Q ) 



i.e. 



Hence, by (3.13), 



St + cSo; > 



C'{U ) 2 

16C^ ' 



So S 1 blowups before t 
Theorem 11.11 



c'(uo)s(o o) wri ich is in 0(1) by (3.14). We complete the proof of 



4. Systems in energy-dependent coordinates 

In this section, by restricting our consideration on smooth solutions, we derive a semi-linear 
system on new coordinates. In the next two sections, we construct the solution for the new 
semi-linear system, then after the reverse transformation we show that this solution is also a 
weak solution of the original system. 
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4.1. Energy-dependent coordinates. We denote 

(4.1) R = (Rx,R 2 ,R 3 ) d = n t + c(ni)n x , S = (Si, S 2 , S3) = f n t - c(ni)n x . 

Without confusion, we still use the letters R and S to denote Riemann variables as in the 
previous section. Then (1.3) can be reformulated as: 

d t Rx - cim^Rx = j^{{c 2 {ni) - l)(\R\ 2 + \S\ 2 ) - 2(3c 2 (m) - i)R ■ S}n x 

+ 2 ^0 j (Ri-S l )R 1 -^(R 1 + S 1 ), 
d t S! + cimfoSt = 44^{(c 2 (n!) - 7 )(| J R| 2 + \S\ 2 ) - 2(3c 2 (nx) - 7 )i? • S}m 
-^(^-SJSi-^Rx + Si) 



(4.2) 



d t R 2 - c(ni)d x R 2 
dtS 2 + c(n 1 )d x S 2 - 
d t R 3 - c{ni)d x R 3 
d t S 3 + c( ni )d x S 3 -- 



2c(m) ' 

= ^{(c 2 ^) - a)(\R\ 2 + \S\ 2 ) - 2(3c 2 (m) - a)R ■ S}n 2 
+^{R2-S 2 )R X - %(R 2 + S 2 ), 

^{(c^i) - a)(\R\ 2 + \S\ 2 ) - 2(3c 2 K) - a)R ■ S}n 2 

-g(fl2-4)Sl-^ + 4), 

= ^{(c 2 ^) -a)(\R\ 2 + \S\ 2 ) - 2(3c 2 (m) - a)R ■ S}n 3 

+^(i?3-S 3 )i?l-^(i?3 + 5 3 ), 

^{(c^i) - a)(\R\ 2 + \S\ 2 ) - 2(3c 2 (m) - a)R ■ S]n 3 



n, 



R-S 
2c(n0 



or lit 



-^(Rz-S^-^R. + S,), 
. R+S 



The energy equation (2.9) equals to 
1 



(4.3) 



^ t {\R\ 2 + \S\ 2 )-\d x [c( ni )(\Rf-\S\ 2 )] 



Pi a 



R + S\ z < 0. 



We define the forward and backward characteristics as follows 

— x ± (s,t, x) = ±c(m(s, x ± (s, t, x))), 



(4.4) 



(Is 
X 



--t = x. 

Then we define the coordinate transformation: 



X 



def 



x~(0,t,x) 



[1 + \R\ 2 (0,y)]dy, and Y 



def 



This implies 
(4.5) 



X t - c(ni)X x 



>x+(0,t,x) 
0, Y t + c{n 1 )Y x = 0. 



[l + \S\ 2 (0,y)]dy. 



Furthermore, for any smooth function /, we get by using (|4.5|) that 
(4.6) 



ft + c( ni )f x = (X t + c( ni )X x )f x = 2c(n 1 )X x f x 
ft - c( ni )f x = (Y t - c(n x )Y x )f Y = -2c( ni )Y x f Y . 



Using (4.6 ), we can transform the directional derivatives dt ± cd x on the (t, x)-coordinates 
into dx and dy on the (X, y)-coordinates, hence get a new semi-linear system on (X, Y)- 
coordinates. We start the calculation from ux and uy . To complete the system, we introduce 
several new variables: 

defl + |i?| 2 defl + |S| 2 



(4.7) 
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and 



(4.8) 



, def 1 
hi = — 

1 + \R\< 



def R 



l + \R\ 2 

def 1 



m = (mi,ni2, m^) 



def 



and 



h 



l + |5| s 



For smooth solutions, \rn\, \h\\ and \h>2\ are all less than 1. However, \p\ and \q\ might 
go to infinity when the possible gradient blowup happens. In the next section, in order to 
estimating them, we will consider the energy equation, where the energy equation can help 
us estimate p and q since X and Y are energy- dependent coordinates. 
We derive the new system in several steps: 



1. Using (4.6) and ( |4.8| ), we have 
(4.9) 







1 



yn 



(d t n - c(ni)d x n) 



2c(m)(-y a 

9jn= . — sv (d t n + c(ni)d x n) 1 



2c{ni 



-m, 



2c( ni )X x 



2c(m) 



2. Then 

d t p-c(n x )d x p = 2{X x )- l [R-(d t R-c(n l )d x R)) 

-(x x )- 2 [d t x x - c( ni )d x X x ](l + \R\ 2 ). 

While using (|4.2l), we have 



R ■ (d t R - c{m)d x R) = - J. J (c 2 (n 1 )-a)(\R\ 2 + \S\ 2 )-2(3c 2 (n 1 )-a)R-S}R-n 

4<r(ni) 



+ ^A(|i?| 2 + \S\ 2 -2R- S)R ini + C ' (m) 



4c 2 (ni) 



— -^R^Rf-R-S) -%(\R\ 2 + R- S), 



while notice that |n| = 1 and d{n{) = ^ 7 e ( ° 1 ' > " 1 ; so that R ■ n = 0, and there holds 

^i?i(|i?| 2 -|5| 2 )-|(|i?| 2 + J R-5), 



(4.10) 



R • (St« - c( ni )d x R) 



which together with (|4.5[) applied gives 



%> - c(ni)d x p 



P 



l + \R\ 2 

from which, and ( 4.6|4.7 ), we infer 

PY= 2c(n!)(-Y^ Pt - c{ni)p 

(4.11) 



ppl(_ Rl{1 + |5|2) + 5l(1 + |£|2)) _ ( |£|2 + ^ 

2c(ni) y 



c (ni, 



i?i 



Si . I-^I 2 -\- R • S 

) -n- 



Similarly, 



2c(m) L2c(m) v l + \R\ 2 l + l^l 2 (l + |S| 2 )(l + |i? 

d t q + c( ni )d xq = 2(-Y x y 1 [S-(d t S + c(n l )d x S)} 

+{YxY 2 (dtY x + c{ ni )d x Y x )(l + \S\ 2 ), 
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then by (4.2), we have 



s ■ (d t s + c( ni )d x s) 



4c 2 (m) 



{(c 2 (ni) - a)(\R\ 2 + \S\ 2 ) - 2(3c 2 (ni) -a)R-S}S-n 



+ 



a — 7 



(\R\ 2 + \S\ 2 -2R- SjSim 



c'(ni) 



2<j ^ )Sl (\S\ 2 -R-S)-^\S\ 2 + R-S), 



4c 2 (m) 

which along with the fact that 5 • n = leads to 

c'(ni) 



(4.12) 



5 • (d t S + c( ni )d x S) = -ppL Sl (\R\i - \S\ 2 ) - ^(\S\ 2 + R-S), 



from which and (4.5), we infer 

q 



d t q + c(m)d x q 



c'( ni ) 

2c(m) 



4c(ni) 



^(1 + 1^2) - Sl (l + \R\ 2 )) - ^(\S\ 2 + R-S) 



Next by (4.6 4.7), we have 



(4.13) 



3. Next, 



pq 



c'im) , R x 



Si \ \S\ 2 + R-S 



2c(m) L2c(n 1 ) v i + |i!| 2 l + \S\ 2 ' (l + |S ? | 2 )(l + |i?| 2 ) 

dth - c( ni )d x h = (I + \R\ 2 )-\d t Ri - cin^dxR!) 

-2(1 + \R\ 2 )- 2 Ri [R ■ (d t R - c{m)d x R) } . 



which together with (4.2) and (4.10) implies 



d t h - c( ni )d x h = (l + ii?! 2 )-^ 



c 2 (ni) -7/|di2 , ,c?i2x 3c2 ( n i)-75 c 



4c 2 (ni) 



Jinx) 



2c 2 (m) 



+(1 + li?! 2 )- 2 ^-^ [ (jRl (i + |5| 2 ) - SxU + \Rf 



+(1 + \R\ 2 )-' 2 R lf i(\R\ 2 + R-S)-(l + li?! 2 )- 1 ^^! + Si). 



Using gggjfr , we obtain 



(4.14) 



+ 



8c 3 (m) 

c'(m) 
Ac 2 (ni) 



[(c 2 (ni) - 7 )(fcx + h 2 - 2/ii/t 2 ) - 2(3c 2 (ni) - 7 )/- m]n x 



- mi) + [24(^2 - fcite + £ • m) - hh 2 - mi/ij . 



4c(m) 



Similarly by (4.8), we have 



a t mi + c( m )^mi =(1 + l^l 2 )- 1 ^! + c(m)d x Si) 

- 2(1 + isi 2 )- 2 ^ [s ■ (d t s + c{m)dj) ] . 
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which together with (4.2) and (4.12) ensures that 



Stmi + cin^mx = (l + ISI 2 )"^ 



lr c 2 (ni)-7 52 32 3c 2 (ni)-7 5 



4c 2 (rai) 



i\R\ z + \s\ 2 ) 



2c 2 (m) 



R-S]m 



- ( ; - \S) 2 r 2 Si^\[(Ri(l + \S\ 2 ) -S 1 (l + \R\ 2 )] 



2c(m) 



+(i + \s\ 2 )- 2 sm\s\ 2 + r.s)-(i + isi 2 )- 1 ^, + s t i. 



from which and (4.6 4.7), we infer 
P 



d x mi 



[(c 2 (ni) - i)(hi + h 2 - 2/ii/i 2 ) - 2(3c 2 (ni) - i)£- m]m 



up 

mip{h - rni) + —, — r [2mi(/n - /ii/i 2 + £ ■ rh) — £\h 2 - mi/ti] . 



8c 3 (m) 

< 4 ' 15 > _ 

4c 2 (ni)""""^ ' 4c(ni) 

Following the same line, we deduce from (4.2) and ( |4.10[ ) that 



d t £ 2 -c(n 1 )d x £ 2 = (l + \R\ 2 )- 1 [ 



lr c 2 (ni)-a 5 2 - 2 3c 2 (ni)-a 5 



4c 2 (m) 



2c 2 (m) 



R-S \n 2 



+(i + \r\ 2 )- 2 r 



c'(m) 
2c(m) 



[R 2 (l + I.SI 2 ) - 5 2 (1 + \R\' 



+(1 + |il| 2 )- 2 i? 2 /x(|i?| 2 + R-S)-(1 + \R\ 2 )- l ^-(R 2 + S 2 ), 



and 



d t £ 3 -c(n 1 )d x £ 3 = (l + \R\ 2 y 1 [ 



lr c 2 (ni)-a 52 - 2 3c 2 (ni)-a 5 



4c 2 (m) 



W + \s\' 



2c 2 (m) 



ii-5 n 3 



+(1 + \r\Y 2 Ri^\ [R S (1 + \S\ 2 ) - S 3 (l + |i?| 2 ) 



2c(m) 



|2\-1/*, 



+(1 + | J R| 2 )- 2 J R 3 /x(| J R| 2 + i? • S) - (1 + l^)- 1 ^^ + S 3 ), 



which together with (|4.6||4.7|) derive that 
dyl 2 



+ 



(4.16) 



+ 



8c 3 (m) 
c'foi) 

4c 2 (m) 
g 

8c 3 (m) 
c'jm) 
Ac 2 (m] 



[(c 2 (ni) - a)(h x + h 2 - 2h x h 2 ) - 2(3c 2 (m) - a)£ • m] n 2 



4c(ni) 



[2^ 2 (/i2 - kite + t- rh) ~ hh 2 - m 2 h l ] . 



[(c 2 (ni) - a) (hi + h 2 - 2h x h 2 ) - 2(3c 2 (ni) - a)£ ■ ml n 3 



£ x q(£ 3 - m 3 ) + [2^3(^2 - /ii/i 2 + f ■ m) - ^2 - m 3 /n] . 



While we deduce from (4.2) and (4.12) that 



d t m 2 + c{m)d x m 2 = (l + ISI 2 )-^ 



2^ir c2 ( n i) - a (|^|2 + i^2n . 3c 2 (ni) -a 



Ac 2 (m) 



2c 2 (m) 



R-S ]n 2 



-(1 + \S\ 2 )- 2 S X ^\[R 2 (1 + \S\ 2 ) - 5 2 (1 + \R\ 2 ) 



2c(m) 



+(1 + |5| 2 )- 2 5 2/ i(|S ? | 2 + i? • S) - (1 + l^l 2 )- 1 !^ + .s',j. 
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and 

d t m 3 + c(m)d x m 3 = (1 + I5I 2 )" 1 f^l'^ m 2 + \S\ 2 ) ~ ^{ni)" * ' ^ 

" (1 + ^~ 2Sl to$ ^ 3(1 + '^ |2) " 53(1 + ^ |2) ] 

+(1 + \s\ 2 y 2 s 3f i(\s\ 2 + r.s)-(i + isi 2 )- 1 ^ + s 3 ), 



which together with (4.6 4.7) implies that 

dxrri2 = —^7 — r T(c 2 (ni) - a) (hi + h 2 — 2h 1 h 2 ) - 2(3c 2 (ni) - a)£ '■ m]n 2 
8c 3 (m) L J 

mip{&2 - m 2 ) + —. — r [2m 2 (hi - hih 2 + £ ■ m) - £ 2 h 2 - m 2 hi] , 



(4.17) 4c ^ ni) 4c(ni) 

9x^3 = g c 3^ n \ [( c2 ( n i) ~~ a )(^i + h 2 - 2hih 2 ) - 2(3c 2 (ni) - a)£ • fh]n 3 

(^l ) ZAP ~~ * 

- 4 c 2( ni ) mi P^ 3 ~~ ms ) + 4 c ( ni ) [ 2m 3(^i - hih 2 +£-m) - £ 3 h 2 - m 3 /ii] . 



4. It follows from (4.10) and (4.8) that 



d t hi - c(m)5 x /ii = - 2(1 + |il| 2 )- 2 [i? • (3 t i? - c(m)^fl) 



(i + l^l 2 )" 2 [^gy^d^l 2 " l^l 2 ) ~K\R\ 2 + R- s)] . 



Then we get by using (4.6 4.7| ) that 

c'K) 



(4.18) 



dyh\ 



4c 2 (nx) " k2) + 2^K) 9/ll(/l2 " ^ 2 + 1 ' ™)' 



Similar calculations together with (4.12) gives 

d t h 2 + c{m)d x h 2 = (1 + |S| 2 )- 2 [^^S^R] 2 - \S\ 2 ) + ^(\S\ 2 + R-S)], 



which together with (4.6 4.7) implies that 
(4.19) 



dxh 2 = 4c2 ^^ pmi(/t 2 - /ii) + 2c( ni ) ph2 ^ hl ~ hlfl2 + ^ ' 
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5. In summary, we obtain 
(4.20) 



dyti = 
d x mi 
b\£ 2 = 
d x m 2 

d Y £ 3 = 

d x m 3 



[(c 2 ( ni ) - 7)(/t! + h 2 - 2h x h 2 ) - 2(3c 2 (ni) - i)t ■ m]m 



c'(m) 



mi) + 



■ [2£ x (h 2 - hih 2 + £-m)- £ x h 2 - m x hi] , 



1 4e 2 (ni) tlC 'V il "" L > 1 4c(ni) 

[(c 2 ( ni ) - 7 )(/ii + h 2 - 2h x h 2 ) - 2(3c 2 (ni) - j)£- m]m 
i^lL™. ^ _ mi ) + ^ [2 mi (fci - hi/12 + m) - tih 2 - mih 

2(3c 2 (ni) -a)£ ■ ml ti-2 
/ii/i 2 + £ • m) - £ 2 ^2 - m^/i-i 



'4?(^) m ^( £ l 



4c(ni) 

[( c 2( ni ) _ a )(/ ll + h2 _ 

+ 4^)^-- 2 ) + 4^7[^( ft 



i h 



[(c 2 (m) - a)(fci + /i 2 - - 2(3c 2 (ni) - a)£ ■ m]n 2 

"' 1 - - m 2 ) + [2m 2 (/ii - /ii/i 2 +/• m) — £ 2 /t 2 - m 2 ^i] . 



[{c 2 { ni ) - a) (Zii + " 2/ii/i 2 ) - 2(3c 2 (ni) - a)£ • m]n 3 



m 3 ) + 



[2£ 3 {h 2 - hih 2 + £-m)- £ 3 h 2 - m 3 h x ] . 



+4c 2 '(m)* 1 ^ v " 3 "Kij t 4c(ni ) 
g^^y [(c 2 K) - a)(fei + /i 2 - 2/ii/i 2 ) - 2(3c 2 (ni) - a)f • m]n 3 

^-mip(* 3 - mg) + jfe [2m 3 (/»i - ^ 2 + * ■ m) - ^ 3 /i 2 - m 3 /ii 



~4c 2 (m)' 



9yn : 
dyhi 

dxh 2 

PY= 2 £,[ 



m, or dx*i 
qhihx 
pm\{h 



2c(m) 
c'(ni) 

4c 2 (m) 
pg r c '( n i) 



2c(m) L 2c(m) 

pq r c'(rai) 
2c(ni) L2c(ni) 



2c(nl) , •' 

2c(n 1 )Qf l i{h 2 - h\h 2 + £■ m), 

(£i - mi) - /i(/t 2 - /ii/i 2 + £ ■ m)] , 
(^i — mi) — n(h\ — h\h 2 + £ ■ m)l . 



4.2. Consistency of variables. Before we prove the global existence of (4.20), we first show 



(1.3) and definitions (4.7) and (4.8) are not used. 



that the various variables introduced for (4.20) are consistent, following from the proposition 



below. Note we only use (4.20) in the proof of this proposition, while the original equation 



Proposition 4.1. For smooth enough data, the following conservative quantities hold: 



(4.21) 



n(X,Y) =m-n(X,Y) = \n(X,Y)\ = l and 



\£(X,Y)\ 2 + h 2 (X,Y) = h 1 (X,Y), \m(X,Y)\' z + h z 2 (X,Y) = h 2 (X,Y) V X, Y, 



2 , l2/ 



as long as 



£ ■ n(X, <p{X)) = m • n(X, ip{X)) = \o.{X,<p(X))\ = l and 
\£(X, <p(X))\ 2 + h 2 (X, <p(X)) = hi(X, tp{X)), 
\m(X, ^{X))\ 2 + h 2 (X, <p(X)) = h 2 (X, <p{X)). 



Here <p represents the initial curve, see (5.1) from the next section. 
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Proof. We first deduce from (4.20) that 
dy [\i\ 2 + h\ - h 



4c 3 (ni) 
(a-7)m 



[(c 2 (ni) - a) (hi + h 2 - 2h 1 h 2 ) - 2(3c 2 (n 1 ) -a)£-rh]£-n 
q[(hi + h 2 - 2h 1 h 2 ) - 2£- fn)£i 



4c 3 (m) 

+ ^T^y^i W? " ^ ™) + 2h i( h i ~ M ~ fa - h 2 )} 
+ 2^)^ - hih 2 + £■ rfi) - £h 2 - mhi] 



+ 



fiq 



2c{m) 

from which and the fact that c/(m) — - — — — 
(4.22) 



(2hj - hi)(h 2 - h x h 2 + £ ■ m), 



c(ni) 



we deduce that 



d Y [\i\ 2 + h\-h l 



4c 3 (m) 



[(c 2 (ni) - a)(hi + h 2 - 2h x h 2 ) - 2(3c 2 (n 1 ) - a)£- m]l ■ n 



+(^K q £i+ m 



•2c 2 {niY 2c(m) 



[h 2 - 2hih 2 + 2£- m]) [\£\ 2 + h\- hi] . 



Similarly it follows from (4.20) that 
dy[£- n] =dy£- n + I- d Y n 



8c 3 (m) 



[(c 2 (ni) - a)(hi + h 2 - 2h x h 2 ) - 2(3c 2 (rai) - a)£- m] (|n| 2 - 1) 



+ J v (hi + h 2 - 2hih 2 ) [(c 2 (ni) - a) + (a - 7 )n 2 ] 
8c? (ni) 1 J 

£ • m[3c 2 (ni) — a + (a — 7 )n 2 ] 



+ 



4c 3 (m) 

Q 

4c 2 (ni 
fiq 



q£l[£ ■ n — m • n] + 



2c(m) 



+ [£(/» 2 - 2/ii/i 2 + 2£ • m) - m/ii] • n, 

which along with the fact that c 2 (ni) = a + (7 — a)n 2 leads to 
(4.23) 

dy[i- n] = - q \(c 2 (ni) - a)(hi + h 2 - 2h x h 2 ) - 2(3c 2 (ni) - a)£- ml (|n| 2 - 1) 
<&cy{ni) 

+ [ ni \ q£i[£- n - m • n] + -r^(h 2 - 2h x h 2 + 2f • m) (f - n) - ^\ (m ■ n). 



4c 2 (m) 



4c(ni) 



4c(ni) 



On the other hand, observe that cVn = 2c (ni) T ^' * s consistent with c?x n = 2e ^^ l.In fact, 
again thanks to (4.20), one has on the one hand 



d x [dyn] 



c'(m) 



8c 3 (m) 

and on the other hand 
c'(ni) 



+ mi)m 



4c 2 (m) 



pq(hi — hih 2 + £ • m)m + 



2c(m) 



8c 3 (m) 



pq(£l + mi)£- —JL—pqfa - h x h 2 + i- m)£+ ^JL^dyt, 



4c 2 (ni) J 



2c(m) 
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which along with the t equations and rh equations of (4.20) shows that (9x[cV n ] = cV[<9x n ]- 
So we can also use the equation dxn = 2c (ni) ^' f rom which and the rh equation of (4.20), we 
get 

dx [rh ■ n] =dxrh ■ n + rh ■ 9x n 
P 



"8c 3 (m) 



[(c 2 (ni) - a){hi + h 2 - 2h 1 h 2 ) - 2(3c 2 (m) - <*)£■ rh] (|n| 2 - 1) 



+ 
+ 

+ 
+ 



P 



[(c 2 (ni) - a)(hi + h 2 - 2h 1 h 2 ) - 2{3c 2 (n 1 ) - a)£ ■ rh] 



8c 3 (m) 

° ~ 7 pn\ [(hi + h 2 - 2h x h 2 ) - 21- 

P 



8c 3 (m) 



m 



4c 2 (m) 



pmi [m • n 



n + 



2c(ni) 



• m, 



MP 



4c(ni) 



[m(/ii - 2/ii/i 2 + 2t ■ rh) - £h 2 ] ■ n, 



which gives rise to 
(4.24) 



P 



• n] = y [(c 2 (m) - a){h x + h 2 - 2h x h 2 ) - 2(3c 2 {n 1 ) - a)£ ■ rh] (|n| 2 - 1) 



(? (n 1 ) > MP ~* 

pmi [£ ■ n — m • n] + — - — -(h\ — 2h\h 2 + 2£ ■ rh)(rh ■ n) 



4c 2 (m) 

While it is easy to observe that 
(4.25) 



4c(ni) 



^[|n| 2 -l] 



pph 2 
4c(ni) 



(£-n). 



c(ni) 



m • n. 



Finally to control the evolution of \fh\ 2 + h 2 — h 2 , we get by applying ( 4.20| ) that 



d x \\rh\ 2 + hi - h 2 



P 



4c 3 (ni) 



[(c 2 (ni) — a)(/ii + h 2 — 2h\h 2 ) — 2(3c 2 (ni) — a)£ ■ rh]fh ■ n 



+ 



Q — 7 

4c 3 (ni) 
c'(m) 



+ /12 — 2h\h 2 ) — 2£ ■ fh]m\n\ 



\2rh(h\ — h\h 2 + £ • rh) — £h 2 — rhhu 
(2h\ - h 2 ){hi - hih 2 + £■ rh) 



MP -> 
H ~, — rwi 

2c(m) 

MP 



2c(m) 



which leads to 
(4.26) 
d x [\rh\ 2 + h 2 2 - h 2 ] 



P 

4c 3 (m) 



[(c 2 (ni) — a)(/ii + /12 — 2h\h 2 ) — 2(3c 2 (n\) — a)£ ■ rh]rh ■ n 



+ 



pmi + ^jF^yfa - 2/l i /l 2 + 2f ■ m]) [\m\ 2 + h 2 - h 2 ] 



■2c 3 (ni 



Summing up (4.22) to (4.26) gives rise to (4.21). This completes the proof of the proposition. 

□ 



singularity and existence to a wave system of nematic liquid crystals 
5. Solutions in the energy coordinates 
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In this section, we prove the existence of the solution for (4.20) with boundary data con- 
verted from (1.5). To avoid the confusion, the rea der sh ould aware that in this section we 
solve variables u, p, q, £, rh, h\ and h% by system (4.20) instead of (1.3) and we do not use 
the definitions (4.7) and (4.8) except when we assign the boundary data. 

The initial line t = in the (t, x)-plane is transformed to a parametric curve 

(5.1) j: Y = <p{X) 

in the (X, Y) plane, where Y = (p(X) if and only if there is an x such that 

X = f*[l + \R\ 2 (0,y)]dy, 

Y = £[l + \3\\0,y)]dy. 

The curve is non-characteristic. We introduce 

def 1 



(5.2) 



(5.3) 



S 



[\R\\0,y) + \S\ 2 (0,y)] dy < 



CO. 



It equals to the number in (1.13) with t = 0. The two functions X = X(x),Y = Y(x) from 
(5.2) are well-defined and absolutely continuous, provided that (1.5) is satisfied. So tp(X) is 
continuous and strictly decreasing on X since X(x) is strictly increasing while Y(x) is strictly 
decreasing. From (5.3) it follows 



(5.4) 

As (t, x) ranges over the domain 
set 



\X + ip(X)\ < 4£ . 

x M, the corresponding variables (X, Y) range over the 



(5.5) 

Along the curve 



{(X,Y); Y>tp(X)}. 



{(X,Y); Y = <p{X)} C 



parametrized by x \-t (X(x), Y(x)j , we can thus assign the boundary data (£, rh, hi, h%,p, q, n) 
G L°° defined by their definition evaluated at the initial data (1.5), i.e., 



n 



n (x), p=l, 



1 



£ = R(Q,x)hi, m = S(0,x)h 2 



(5.6) 



hi 



1 



1 + |i2| 2 (0,x) 
1 

l + |5P(0,x)' 



and 



where 



R(0, x) = rii(x) + c(ni (x))rio(x), 5(0, x) = m(sc) - c(ni (x))ri (x) . 
We consider solutions to the boundary value problem ( |4.20 )(5.6)( |1.4| ). 



Theorem 5.1. The problem (4-20)(5.6)(1.5)(1.4) has a unique global solution defined for all 
(X,Y)en+. 
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Sketch of Proof. Noticing that all equations in (4.20) have a locally Lipschitz continuous 
right hand side, the construction of a local solution as fixed point of a suitable integral 
transformation is straightforward. Note the consistency condition (9x[<9y n ] = &Y [dx n ] proved 
in the previous section shows that we can use either the equation for dyn or the one for dxn 
to find the solution. 

To make sure that this solution is actually defined on one must establish a priori 
bounds, showing that the solution remains bounded on bounded subset of Q + . 
By ( |4.21[ ), we have 

(5.7) 



Ml-/ii)>0, h 2 (l-h 2 )>0. 



Thus hi,h 2 are bounded between zero and one, and \£\ and \fh\ are both uniformly bounded 
by ( pH >. 

By the p and q equations in (4.20) and (4.21), we have 
(5.8) 



py + qx 



mh\h 2 + fh 2 < 
2c(n\) hi h 2 ~ 



which implies that 

r-X 



[ p{X',Y)dX' + [ q(X,Y')dY' < X - tp~ l (Y) + Y - <p(X) 
J v -i(Y) J V {X) 

where </? _1 denotes the inverse of ip, following an integration over the characteristic triangle 
with vertex (X,Y). Thus, by the energy assumption (5.3), we find 



(5.9) 



f X ,Y 

/ p{X',Y)dX'+ q{X,Y')dY' <2(\X\ + \Y\+4£ ). 

J v -i(Y) Jip{X) 



Integrating the p equation in (4.20) vertically and use the bound on q from (5.9), we find 
p(X,Y) 



exp 



(5.10) 



«x) q Wiy [~ C M^ ~ m ^ ~ M ~ h ^ + 1 -)] dY '} 

< ex P {c f* x) q(X,Y>)dY>} 

< exp{2C (|X| + |y|+4£:o)}. 
Here Co represents a finite number. Similarly, we have 

(5.11) q(X, Y) < exp{2C (|X| + \Y\ + 45 )}- 

Relying on the local bounds (5.10) ( [57TT] ) , the local solution to (4.20)(5.6)(1.5)(1.4) can be 
extended to Sl + One may consult paper [3J for details. This completes the sketch of the proof. 

Corollary 5.1. If the initial data (no,ni) are smooth, the solution U := (n,p, q, £, fh, hi, h 2 ) 
of (4-20)( 5.6)(1.5)(1.4) is a smooth function of the variables (X, Y). Moreover, assume that 
a sequence of smooth functions (n^n^)^! satisfies 



no 



n 



Qjx 



(no) a 



n. 



ni 



uniformly on compact subsets of If 
solutions: 

(<^g\^m 4 ,/4,/4) 
uniformly on bounded subsets offl + . 



Then one has the convergence of the corresponding 
U 



SINGULARITY AND EXISTENCE TO A WAVE SYSTEM OF NEMATIC LIQUID CRYSTALS 



19 



6. Inverse transformation 
By expressing the solution n(X, Y) in terms of the original variables (t, x), we shall recover a 



solution of the Cauchy problem p.3|~(1.4). This will prove Theorem 1.2 except the estimate 



£(t) < £q which will be proved in the next section. 



Using (4.6) by letting / = t or x, we have equations: 



/ C1l . Phi , qh 2 -qh 2 

(6-1) tx = -z— , ty = — — , X X = —, Xy 



2c 2c 2 ' 2 

In fact, we only need one of the two equations of tx and ty, which are consistent since 



txY = tyx- The same is true for x. We integrate (6.1) with data t = 0, x = x on 7 to find 
t = t(X, Y),x = x(X, Y), which exist for all (X, Y) in Q+. 

We need the inverse functions X = X(t,x),Y = Y(t,x). The inverse functions do not 
exist as a one-to-one correspondence between (t, x) in R + x R and (X, 1") in Q + . There may 
be a nontrivial set of points in Q + that maps to a single point (t,x). To investigate it, we 
find the partial derivatives of the inverse mapping, valid at points where hi 7^ 0, h 2 7^ 0, 

(6-2) X t = ^T> y * = ^T' X - = "^' Y * = ~\- 

phi qh 2 phi qh 2 



Thus (4.5) holds and so does (4.6) for our solution. 

We first examine the regularity of the solution constructed in the previous Section. Since 
the initial data (no) x ,ni etc. are only assumed to be in L 2 , the functions U may well 
be discontinuous. More precisely, on bounded subsets of the + , the solutions satisfy the 
following: 

• The functions £i,i 2 , ,£3, hi,p are Lipschitz continuous w.r.t. Y, measurable w.r.t. X. 

• The functions mi,m 2 ,m,3, h 2 , q are Lipschitz continuous w.r.t. X, measurable w.r.t. Y. 

• The vector field n is Lipschitz continuous w.r.t. both X and Y. 

In order to define n as a vector field of the original variables t,x, we should formally 
invert the map (X, Y) h-> (t, x) and write n(t,x) = n(X(t,x) , Y(t, x)). The fact that the 
above map may not be one-to-one does not cause any real difficulty. Indeed, given (t*,x*), 
we can choose an arbitrary point (X*,Y*) such that t(X*,Y*) = t* , x(X*,Y*) = x* , and 
define n(t*,x*) = n(X* ,Y*). To prove that the values of n do not depend on the choice 
of (X*,Y*), we proceed as follows. Assume that there are two distinct points such that 
t(X 1 ,Y 1 ) = t(X 2 ,Y 2 ) = t*, x(X u Y-y) = x(X 2 , Y 2 ) =x*. We consider two cases: 
Case 1: X\ < X 2 , Y% < Y 2 . Consider the set 

I> := {(X,Y); x(X,Y) < x*} 



and call dT x * its boundary. By ( |6.1| ), x is increasing with X and decreasing with Y . Hence, 
this boundary can be represented as the graph of a Lipschitz continuous function: X — Y = 
(j)(X + Y). We now construct the Lipschitz continuous curve 7 consisting of 

• a horizontal segment joining (X\,Y\) with a point A = (Xa,Ya) on dT x *, with 
Y A = Y X , 

• a portion of the boundary dT x * , 

• a vertical segment joining {X 2 ,Y 2 ) to a point B = (Xb,Yb) on dT x *, with Xb = X 2 . 



Observe that the map (X, Y) 1— > (t, x) is constant along 7. By (6.1 ) this implies hi = on the 
horizontal segment, h 2 = on the vertical segment, and hi = h 2 = on the portion of the 
boundary dT x * . When either hi = or h 2 = or both, we have from the conserved quantities 



(4.22) that £ = or m = or both, correspondingly. Upon examining the derivatives of n in 
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(4.20), we have the same pattern of vanishing property. Thus, along the path from A to B, 
the values of the components of n remain constant, proving our claim. 
Case 2: X\ < X2, Y\ > I2. In this case, we consider the set 



IV := {(X,Y); t(X,Y) <f } 



and construct a curve 7 connecting (X\, Y{) with (X2,Y2) similarly as in case 1. Details are 
entirely similar to Case 1. 

We now prove that the function n(t,x) = n(X(i, x), Y(t, x)) thus obtained are Holder 
continuous on bounded sets. Toward this goal, consider any characteristic curve, say t 1— > 
x + (t), with dx + /dt = c(rai). By construction, this is parametrized by the function X 1— > 
(t(X , Y), x(X, y)) , for some fixed Y. Using the chain rule and the inverse mapping formulas 
d6.2|), we obtain 



n t + cn x = n x (X t + cX x ) + n Y (Y t + cY x ) = 2cX x n x . 



Thus we have 



Jo K + 



cn. 



dt 



J^(2cX x \n x \f(2X t )- 1 dX 



(6.3) 



1 pW \2phi 
ph x 2c > 2c 



= Ix V^dX < f£i- c dX < C T , 

for some constant C T depending only on r. Notice we have used \£\ 2 < hi, which follows from 
( |4.21 ). Similarly, integrating along any backward characteristics 1 1— > x~(t) we obtain 



dX 



(6.4) 



cnJ dt < C T . 



Since the speed of characteristics is ±c(ni), and c(ni) is uniformly positive and bounded, the 
bounds ( 6.3 )-( 6.4 ) imply that the function n = n(i, x) is Holder continuous with exponent 1 /2. 
In turn, this impli es th a t all characteristic curves are C 1 with Holder continuous derivative. 
In addition, from (6.3)-(6.4) it follows that R, S at (4.1) are square integrable on bounded 
subsets of the t-x plane. However, we should check the consistency that R, S at (4.1) are 
indeed the same as recovered from (4.8). Let us check only one of them, R = l/h\. We find 



n t + cn x = 2cX x n x = 2c 



1 P £ _ e 

phi 2c hi 



Finally, we prove that n satisfies the equations of system (1.3) in distributional sense, 
according to (iii) of Definition 1.1. We note that 

r c 2 n x ] + 2[i(j)tndxdt 



;/2 



(6.5) 



// 4> t [(n t + cn x ) + (n t - cn x )] 
-ccp x [(nt + cn x ) - (n t - cn x )] 



2fj,(ptn dxdt 



J J [fa - C(p x ] (n t + cn x ) 

+ // [& + c4>x] (n* - cn x ) + 2^4> t u.dxdt 

Jf [4>t — c4>x] R dxdt + J J [(fit + c(fi x ] S dxdt + J J 2/^n dxdt . 
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By ( |4.6[ ), this is equal to 

(6.6) J J [ - 2cY x cp Y R + 2cX x <p x S + 2/ic(X E x - Y x <t>y)n] dxdt . 
Using the Jacobian 
(6.7) 
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d(x,t) -pqhxhi 



d(X,Y) 2c 



derived by (6.1), and the inverse (6.2), it is equal to 

——K(pY + ——bcpx + —ftPx + —r-PY) 
qh 2 phi phi qh 2 



pqhih 2 

2c 



dXdY 



(6i 



J J [phiR4>Y + qh2S(j)x + ^qh2^4>x + Atp/iin0y] dXdY 
Y + qfh(f)x + fJ,qh2ii(f)x + /up/iin^y] dXdY 
[ ~ (pt)Y - (qrn) x - n(qh 2 n) x - n(phin) Y ](j> dXdY. 



Thanks to (4.1), (4.8) and (4.20), we have that the first component of the integrand equals 



to 



4>[-{pii)Y ~ (qmi) x ~ n(qh 2 n)x - n(phin) Y 

, pq 

'4c 3 



</>f| [(c 2 - 7 )(/»i + h 2 - 2h 1 h 2 ) - 2(3c 2 - j)£- m] n x 



pqhih 2 r c 2 - 7 1 1 
L pqhih 2 r c 2 - 7 - 



3c 2 — 7 ^ • m-. 



1"2 



2c 



3c 2 -7 



R-S] 



L pqhih 2[ c 2 -7 2 2 2n 3c 2 -7 
-2— (|n t | +c In^l ) 



2<j) 



2c 

pqhih 2 , , l2 | 2 _ 2 , 



c 2 |n 



2c 



|n t | 2 + (2c 2 - 7 )|n x | 2 m, 



second and third equations of system (1.3) in distributional sense. 



which implies that the first equation in (1.3) holds in integral form. Similarly, n satisfies the 



7. Upper bound on energy 



We convert the energy equation (2.9) formally to the (X, Y)-plane to look for the upper 
bound of energy. 



Recall that the energy equation (2.9) can be written as (4.3) in terms of \R\ and \S\. By 



the variables (|4.8|), we rewrite the equation (|4.3|) as 

1 



(7.1) 



We can write the 1-form 



+ 



1 




[V - 






4/i 2 


-0r 


4% 




X 



< 0. 
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as 



(7.3) 

by the formula 
(7.4) 



dX 



/(I - h 2 ) 



dt 
dx 



t x dX + tydY 
xxdX + xydY 



By (4.20) then by (4.21), we have 



1 



[{p(l-h 1 )) Y +[q(l-h 2 )) 



Xi 



(7.5) 



dY, 



ph 



^ dX + ^dY 



^dX 



<^dY. 



pqp 

' 2c 
pqp 

' 2c 



(hi + h 2 - 2hih 2 + 2£-m) 
(\£ + m\ 2 + (hi-h 2 ) 2 ) 



< 0. 




Figure 2. Energy conservation 



We use ( 7.3 ) to prove that the 

(7.6) S(t) < £ for any t > 0. 
Fix r > 0, and r » 1. Define the set 

(7.7) n+ := {(X,Y); < t(X,Y) < r, X < r , Y < r\. 

See Figure^ where segment A t Aq is where Y = r while segment B t Bq is where X = r. By 
construction, the map (X, Y) \— > (t, x) will act as follows: 



A T h-> (r, a) , 
: c. : 



B r ^(T,b) 



B ^ (0, c) , >-> (0, d) , 



for some a < b and d < c. By the divergence theorem and (7.2), the integral of the form (7.3) 
along the curve A T — > B T — > Bq — > Aq is less than or equal to zero. Integrating the 1-form 
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(7.3) along the boundary of 0+ we obtain 



(7.8) 





4 


9(1 


-h a ) 




4 


9(1 


-h 2 ) 



dY ~ Jaoa/"^ dX - J BoBt 
dY 



= til\\n t \ 2 {0,x) + c 2 (n 1 {0,x)) |n x .| 2 (0,x) 
On the other hand, we use (|7.4|) to compute 



dx . 



(7.9) 



|n t | 2 (r,x) + c 2 (ni(T,x)) \ii x \ 2 (t,x) 
p{\ - hi 



dx 



-dX 



q(l-h 2 ) 



dY < 8 . 



A T B T n{h 2 ^0} 



Letting r — > +oo in (7.7), one has a — > — oo, b — > +oo. Therefore (7.8) and (7.9) together 
imply (7.6). 

8. Regularity of trajectories 

8.1. Lipschitz continuity. In this part, we first prov e the Lipschitz continuity of the map 
1 1—?- (rti, ri2, n^){t, •) in the L 2 distance, stated in (1.12). For any h > 0, we have 



n(i + h, x) — n(i, x) = hi tit(t + rh, x) dr. 

Jo 



Thus 

(8.1) \\rii(t + h,x) — rii(t,x)\\ L 2 < h 



/ limit 

Jo 



+ Th,-)\\ L 2dT <hyj2S , i = l~3. 



8.2. Continuity of derivatives. We prove the continuity of functions t h-> (rii t , n 2 t, n^tjit, •) 
and t i y (rii x ,ni x ,ni x )(t, ■), as functions with values in LP, 1 < p < 2. (To keep tradition, 
we use the exponent p here at the expense of repeating one of our primary variables.) This 



will complete the proof of Theorem 1.2 



We first consider the case where the initial data (no) x and ni are smooth with compact 
support. In this case, the solution n = n(X,Y) remains smooth in f2 + . Fix a time r > 0. 
We claim that 



(8.2) 

where 

(8.3) 



d_ 
df 



nit,-) 



n*(r,-) 



t=T 



pi c qm c £ rh 

m(r,s) ■■=^ x ^^ = YcWi + ^W2 = ^hi + 2-h 2 



Notice that (8.3) defines the values of nt(r, •) at almost every point By the inequality 

d7.6|), we obtain 



(8.4) 



/ \n t (T,x)\ 2 dx < 2£{t) < 2£ . 

JR 
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To prove (8.2), we consider the set 

(8.5) T T :={(X,Y)\t(X,Y)<r}, 

and let jr be its boundary. Let e > be given. There exist finitely many disjoint intervals 
[di, bi] C R, i = 1, . . . , N, with the following property. Call Ai,Bi the points on 7 r such that 
x{Ai) = a%, x(Bi) = h. Then one has 

(8.6) mm {h 1 (P),h 2 (P)} < 2e 

at every point P on 7 T contained in one of the arcs AJBi, while 

(8.7) h 1 (P)>e, h 2 (P)>e, 

for every point P along 7 T! not contained in any of the arcs AiBi. Call J := U\<i<jsr[a,i, bi], 
J' = M\J, and notice that, as a function of the original variables, n = n(t, x) is smooth in a 
neighborhood of the set {r} x J' . Using Minkowski's inequality and the differentiability of n 
on J', we can write, for i = 1 ~ 3, 

lim/j^o ji IJ R ni(r + h,x) - ni(r,x) - hn it (T,x) dx) 



< lim 



h->0 h \ J J 



nj(r + h,x) - rii(T,x) 



v \ 1 /p 
dx 



+ {fj\ n it(i~,x)\ dx) 
We now provide an estimate on the measure of the "bad" set J: 
meas(J) ./.,"- /.,./.., , 

2 >:../; 



Ijd- = Ej AlB ,^dX-^dY 



< 



l-2e 



AiBt 2 UA 2 al 



<• 4e r p(l-fei) j Y 9(1-^2) jy <■ 4e c 
^ T=2i J Tt . — 4 — dA 4 — aYSj^to- 

Notice that dt = on y T , so the two parts of the integral are actually equal. Now choose 
q = 2/(2 — p) so that | + - = 1. Using Holder's inequality with conjugate exponents 2/p and 
q, and recalling (8.1 ), we obtain for any i = 1 ~ 3 

p 

J j riiij + h,x) — ni(r,x) dx 

< meas ( J) 1 ^ ■ (fj nj(r + h, x) — nj(r, x) 

< meas(J) 1/g • M|rai(r + h,-)- n^r, OHxa) 

< meas(J) 1 / 9 • (/i 2 [2f ] ^ 



2 \P/ 2 

dx 



\P/2 



Therefore, 



(8.10) 



lim sup — 

/i->0 « 

4e 



rij(r + h,x) - th(t, x) 



dx 



i/p 



1 1/2 
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In a similar way we estimate 

(r, x)\ p dx < [meas (J)] 



n 



a 


nu(T,x) 







2 \P/ 2 

dx I 



Ul) 



c \ i/p 

y |ni t (r,x)| p da;J < meas ( J) 1 ^ 9 • [2£ ] p/2 . 



Since e > is arbitrary, from ( |8.8[ ), (8.10) and ( 8.11| ) we conclude 



5.121 



lim — 

h^o h 



riiir + h,x) — rij(r, x) — /inj t (r, x) 



dx 



i/p 



0. 



The proof of continuity of the map t i— > nj t is similar. Fix e > 0. Consider the intervals 
[aj, 6j] as before. Since n is smooth on a neighborhood of {r} x J', it suffices to estimate 

limsup^o J |ni t (r + /i, x) — ni t (r,x)\ p dx 



< limsup^o jj |n if (r + - nj t (r,x)| p dx 



< limsup^Q [meas (J)] • ^Jj ni t (r + h, x) —rii t (T,x) dxj 

< limsup^o [j^£q 



\ P/2 



I 1/9 



|n it (r + /i, + ll n ^( r '")|| L 2 



Since e > is arbitrary, this proves continuity. 

To extend the result to general initial data, such that (no)^, ni = n t \t=o € L 2 , we use 
Corollary 5.1 and consider a sequence of smooth initial data, with (ng)^,]!^ G C^°, with 
rig — > no uniformly, (ng)a; — > (iiq) x almost everywhere and in L 2 , — > ni almost everywhere 
and in L 2 . 

The continuity of the function 1 1— > n x (t, •) as maps with values in L p , 1 < p < 2, is proved 
in an entirely similar way. 
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